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Ground state of the periodic Anderson model on a triangular lattice is systematically in- 
vestigated by the mean- field approximation. We found that the model exhibits two different 
types of partially disordered states: one is at half filling and the other is at other commensu- 
rate fillings. In the latter case, the kinetic energy is lowered by forming an extensive network 
involving both magnetic and nonmagnetic sites, in sharp contrast to the former case in which 
the nonmagnetic sites are rather isolated. This spatially extended nature of nonmagnetic sites 
yields a metallic partially-disordered state by hole doping. We discuss the mechanism of the 
metal-insulator transition by the change of electronic structure. 
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Geometrical frustration is the conflict of competing in- 
teractions originating in the structure of underlying lat- 
tice. 1 ' 2 The competition leads to degeneracy in the low 
energy states, which is the source of fascinating phenom- 
ena, such as unusual orderings and exotic ground states. 
One of the peculiar orderings under geometrical frustra- 
tion is a partial disorder (PD), which is characterized 
by coexistence of magnetic order and paramagnetic or 
nonmagnetic sites. This offers an intriguing example of 
self-organization to relieve geometrical frustration. 

PD was first theoretically proposed in the localized 
Ising spin model on a triangular lattice, 3 and later ex- 
perimentally observed in ^-electron compounds. 4, 5 In the 
localized spin systems, PD is mainly driven by the en- 
tropic effect, and hence it is stable only at finite temper- 
atures. On the other hand, a metallic PD was also found 
in some / and d-electron compounds. 6 8 In stark con- 
trast to the insulating PD, the metallic ones persist down 
to the lowest temperature, which suggests that itinerant 
electrons play a key role in stabilizing the PD states. 
In particular, in the /-electron compounds where itiner- 
ant electrons are coupled with localized / moments, the 
Kondo singlet formation 9, 10 is expected to be relevant 
to the vanishment of moments at the nonmagnetic sites. 
In spite of extensive studies, the stabilization mechanism 
and the nature of these metallic PD states are not fully 
understood yet. 

In the present study, we theoretically examine the pos- 
sibility of PD in itinerant electron systems, with focusing 
on the role of coupling between itinerant and localized 
electrons. Similar attempts were made for an effective 
pseudospin model, 11,12 the Kondo necklace model, 13 15 
the Kondo lattice model, 14, 16 and the periodic Anderson 
model 17, 18 on frustrated lattice structures. However, the 
obtained PD states have been all insulating so far, al- 
though metallic PD states are experimentally observed. 
It is highly desired to establish a theory to describe 
metallic PD states not only for understanding the ex- 
perimental results but also for paving the way to explore 



unusual magnetotransport properties associated with the 
peculiar coexistence of magnetic and nonmagnetic sites. 

In this Letter, we systematically investigate the ground 
state of the periodic Anderson model on a triangular lat- 
tice. We identify new commensurate PD states, which 
are distinct from the one previously found at half fill- 
ing. 17 The striking feature is that the nonmagnetic sites 
contribute to the kinetic energy gain by participating in 
an extensive network with the magnetic sites, whereas 
they are rather isolated in the half- filling case. Reflect- 
ing the peculiar feature, we successfully obtain a metallic 
PD state by hole doping. We discuss the nature of the 
metal-insulator transition in detail. 

We consider the periodic Anderson model on a trian- 
gular lattice, whose Hamiltonian is given by 



n 



= ~t E (i c F + H^-^(cU + Ec.) 



^ 1 fia fio- • 



(i) 



The first term represents the kinetic energy of conduc- 
tion electrons, the second term the onsite c-f hybridiza- 
tion between conduction c and localized / electrons, the 
third term the onsite Coulomb interaction for localized 
electrons, and the fourth term the atomic energy of local- 
ized electrons. The sum of is taken over the nearest- 
neighbor sites on the triangular lattice. Hereafter, we 
take t = 1 as an energy unit. 

In order to study the ground state of model Eq. (1), 
we employ the standard Hartree-Fock approximation for 
the Coulomb U term, which preserves the SU(2) sym- 
metry of the interaction term. 17, 18 In the calculations, 
we adopt a three-site unit cell to determine the global 
phase diagram. We also check the stability of the relevant 
phases by using a twelve-site unit cell, which accommo- 
dates longer-period orderings, such as period four (2 x 2) 
and six (2 x 3). 

We performed the calculations mainly at the com- 
mensurate fillings, n = 2m/3 with an integer m (n = 
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Fig. 1. (color online). Ground-state phase diagram of the peri- 
odic Anderson model at (a) 2/3 filling and (b) half filling. Eq 
is taken at —4 and —1, respectively, (b) is partly taken from 
ref. 17 for comparison. Schematic picture of the ordering pat- 
tern is shown for each phase. The sizes of the circles reflect local 
electron densities, and the arrows represent local spin moments. 
PDI, AFM, NM, KI, and FM stand for the PD insulator, AF 
metal, nonmagnetic metal, Kondo insulator, and ferromagnetic 
metal, respectively, (c)-(e) V dependences of (c) the magnitude 
of local spin moment, (d) local charge density, and (e) energy gap 
at 2/3 filling for U = 2 and Eq = —4. In (c) and (d), the data are 
plotted for c and / electrons at each sublattice separately, and 
the dot-dashed and dashed curves represent the mean values of c 
and / contributions, respectively. The solid lines (symbols) show 
the results of the calculations with the three-site (twelve-site) 
unit cell in (c)-(e). 



^2ia( c ia c io- + fiafioj/N, where N is the total number of 
sites), and the doped region around them. As a result, in 
addition to the PD state previously found at half filling 
(n = 2), 17 we obtained new PD states at 2/3 and 8/3 
fillings. In the following, we focus on the result at 2/3 
filling as the nature of the PD state is similar to that at 
8/3 filling. 

Figure 1(a) shows the ground-state phase diagram at 
2/3 filling. The phase diagram is divided into two regions: 
magnetically ordered region for U > V and nonmagnetic 
region for U < V. The magnetic and charge states as 
well as the energy gap are shown in Figs. 1(c)- (e). 

In the region for U < V, the system is in a nonmag- 
netic metallic state. The magnetic moments vanish at all 
the sites in both c and / components [Fig. 1(c)], and the 
system is a gapless metal [Fig. 1(e)]. On the other hand, 
a 120° antiferromagnetic (AF) metallic phase appears 
in the narrow window next to the nonmagnetic metallic 
phase. In this phase, the c and / moments constitute a 
three-sublattice noncollinear 120° order, as schematically 
shown in Fig. 1(a). 

While increasing U or decreasing V, the system 
changes from metal to insulator with exhibiting a PD, 
as shown in Fig. 1(a). 19 In this phase, one of the three 



Fig. 2. (color online). V dependences of (a),(d) the kinetic energy 
for each bond and (b),(e) hybridization energy at each site. See 
the definitions in the text. The data are taken at (a),(b) U = 2, 
E = -4, n = 2/3 and (d),(e) U = 2, E = -1, n = 2. The 
symbols are common to those in Figs. l(c)-(e). (c),(f) Schematic 
pictures of the corresponding PD states; the thick lines show the 
bonds on which the kinetic energy gain is larger. 



sublattices becomes nonmagnetic, while the remaining 
two sublattices retain magnetic moments, which show a 
collinear AF order [Fig. 1(c)]. The moments are larger 
for localized electrons than for conduction ones, and the 
magnitude of each moment takes the same value be- 
tween the two sublattices. The PD is accompanied by 
charge disproportionation, as shown in Fig. 1(d); the lo- 
cal charge density becomes higher at the magnetic sites 
than the nonmagnetic sites. The transition from the 120° 
AF metal to the PD insulator is of first order, while the 
charge gap grows continuously, as shown in Figs. 1(c)- (e). 
We obtain essentially the same results by the calculations 
for the twelve-site unit cell, as plotted in Figs. 1(c)- (e). 

A similar three-sublattice PD insulating state was ob- 
tained by the authors at half filling between the 120° AF 
metallic phase and nonmagnetic Kondo insulating phase, 
as shown in Fig. 1(b); the phase diagram is taken from 
ref. 17 for comparison. However, there are several con- 
trastive aspects between the two PD phases. First of all, 
the two PD states have different structures in the spin 
and charge ordering patterns. As shown in the schematic 
pictures in Figs. 1(a) and 1(b), the c and / moments at 
each magnetic site are aligned parallel to each other in 
the 2/3-filling case, while antiparallel in the half- filling 
case. In addition, the local charge density is higher at 
the magnetic sites than the nonmagnetic site in the 2/3- 
filling case, whereas the tendency is opposite in the half- 
filling case. 

More interestingly, the two PD states show the distinct 
nature of spatial coherency. This is clearly observed by 
calculating the kinetic energy for nearest-neighbor bonds 
(ij), = — t(%2 a c la c jo- +H.c), and the hybridization 
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Fig. 3. (color online). Ground-state phase diagram for carrier 
doping near (a) 2/3 filling and (b) half filling. The data are taken 
at (a) U = 2, E = -4 and (b) U = 2, E = -1. (b) is taken 
from ref. 18 for comparison. Schematic picture of the ordering 
pattern is shown for each phase. PDM and AFI stand for the 
PD metal and AF insulator. PS represents the phase separated 
region. In (a), the PDI and 120° AFI states exist only at n = 2/3 
and n = 1, respectively. In (b), the PDI and KI states exist only 
at n = 2. 

energy at site z, E l v = — V(^2 a c\ G f%a + H.c). As shown 
in Fig. 2(d), at half filling, the kinetic energy gain at the 
bond between a magnetic site and its neighboring non- 
magnetic site is smaller than that between the magnetic 
sites in the PD state. This means that the nonmagnetic 
sites are weakly connected with the magnetic sites which 
form a honeycomb subnetwork of large kinetic energy 
gain, as schematically shown in Fig. 2(f). Corresponding 
to this real-space separation, the c-f hybridization gains 
a larger energy at the nonmagnetic sites than the mag- 
netic sites, as shown in Fig. 2(e). These indicate that the 
c-f hybridization at the nonmagnetic site plays a crucial 
role in stabilizing PD at half filling. 

In contrast, in the PD state at 2/3 filling, the nonmag- 
netic sites are not separated from the magnetic sites, but 
rather form a coherent network with magnetic sites, as 
schematically shown in Fig. 2(c). This is clearly seen in 
the results plotted in Figs. 2(a) and 2(b); the kinetic en- 
ergy gain becomes larger for the bonds including non- 
magnetic sites, and the hybridization energy gain be- 
comes smaller for the nonmagnetic sites compared to the 
magnetic sites. In addition, the hybridization energy gain 
is overall smaller than that at half filling. These suggest 
less importance of the local c-f hybridization for PD at 
2/3 filling compared to half filling. This is reasonable 
since the average /-electron density at the nonmagnetic 
site is typically around 0.5, which is too small to gain the 
energy by forming a c-f bonding state. 

The results above indicate that the nature of spa- 
tial coherency is clearly different between the PD states 
at half filling and 2/3 filling; the nonmagnetic sites are 
rather spatially isolated in the former, whereas they form 
an extended network with magnetic sites in the latter. 
This implies that, when going beyond the mean-field 
level, the nonmagnetic sites form spatially extended sin- 
glet states in the 2/3-filling PD state, whereas they form 
local Kondo singlets in the half-filling case. 

Reflecting the different coherency, the two PD insula- 
tors exhibit distinct responses to carrier doping. Figure 3 
shows the results of the ground-state phase diagram for 



both hole and electron doping to 2/3 filling and half fill- 
ing. The phase diagrams are obtained by using the grand 
canonical ensemble method. 18 

In the 2/3-filling case, hole doping to the PD insula- 
tor induces a metal-insulator transition to a PD metallic 
state, as shown in Fig. 3(a). The PD metallic state be- 
comes stable in a wider range of doping while increasing 
V, whereas it disappears at V > 2.25. On the other hand, 
the electron doping immediately leads to a phase separa- 
tion to the 120° AF metallic state. The situation is quite 
different in the half-filline; case; as shown in Fig. 3(b) 
(taken from ref. 18 for comparison), both hole and elec- 
tron dopings lead to a phase separation and no PD metal- 
lic state is obtained. 

To understand the distinct responses to carrier dop- 
ing between the half-filling and 2/3-filling cases, let us 
discuss the nature of the PD state and metal-insulator 
transition from the viewpoint of the electronic structure. 
Figure 4(a) plots a typical local density of states for the 
PD insulating state at 2/3 filling. The results show the 
contributions from the nonmagnetic (upper half in the 
figures) and magnetic sites (lower half) separately. There 
is a small energy gap at the Fermi level, as shown in the 
inset of Fig. 4(a). By hole doping, the gap closes as shown 
in Fig. 4(b), while the PD magnetic structure persists 
even in the metallic state. 

The robustness of the PD metallic state near 2/3 fill- 
ing can be attributed to the electronic structure near the 
Fermi level. In the 2/3-filling case, the density of states 
just below the Fermi level dominant ly consists of the con- 
tribution from magnetic sites [inset of Fig. 4(a)]. Hence, 
doped holes enter mainly into the magnetic sites. They 
reduce the ordered magnetic moments, but do not dis- 
turb the c-f hybridized states at the nonmagnetic sites 
seriously; consequently, the PD state persists as a metal- 
lic state upon hole doping. On the other hand, for elec- 
tron doping, carriers are mainly doped into the nonmag- 
netic sites as the density of states just above the Fermi 
level is dominated by the nonmagnetic contribution [in- 
set of Fig. 4(a)]. In this case, the doped electrons directly 
disturb the c-f hybridization, and hence, the PD state 
is unstable toward the phase separation to magnetically- 
ordered state. 

In contrast, in the half-filling case, the density of states 
both just below and just above the gap is dominated by 
the contribution from nonmagnetic sites, i.e., the gap 
has a strong c-f hybridization nature. 17 In this case, 
both hole and electron dopings affect the nonmagnetic 
sites and reduce the energy gain by the c-f hybridiza- 
tion. This gives rise to moment formation at the nonmag- 
netic sites, leading to the instability of PD state toward 
magnetically-ordered states. 

The above picture is consistent with the contrastive 
spatial coherency found in Fig. 2. At half filling, doped 
carriers concentrate on the nonmagnetic sites which are 
rather separated from the magnetic subnetwork. It is 
not easy to make the system metallic by the doped c- 
f hybridized states localized at the nonmagnetic sites; 
it seems necessary to reconstruct the magnetic order by 
inducing magnetic moments at the nonmagnetic sites, 
which immediately destroys PD. In contrast, the c-f hy- 
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filling, where the nonmagnetic sites are rather isolated 
from the magnetic honeycomb subnetwork. Reflecting 
the characteristic coherent nature, a metallic partially- 
disordered state is realized by carrier doping. The mecha- 
nism of the metal-insulator transition was analyzed from 
the viewpoint of the electronic structure. Our results will 
offer a good starting point to pursue the magnetotrans- 
port specific to the peculiar coexistence of magnetic and 
nonmagnetic sites. 
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Fig. 4. (color online). Local density of states of conduction and 
localized electrons for (a) PDI at n = 2/3 and (b) PDM at 
n ~ 0.617. The data are taken at U = 2, V = 2, and E = -4. 
In the figures, the upper (lower) half shows the result for the 
nonmagnetic (magnetic) sublattice site(s) 1 (2 and 3). The data 
show the electronic structures near the Fermi level indicated by 
the thin vertical lines; the bands extend in the energy range 
approximately from —7 to 4. In (a), the inset plots the enlarged 
figure in the vicinity of the Fermi level. 



bridized state is spatially extended including the mag- 
netic sites at 2/3 filling. Carriers are not fatal to PD 
state, if they are mainly doped into the magnetic sites. 
In this case, a metallic conductivity is easily acquired 
by the carriers doped on the extended coherent network. 
We expect that the PD metallic state remains stable even 
when we include quantum fluctuations beyond the mean- 
field approximation. This is because fluctuations may 
enhance the singlet formation by the c-f hybridization 
around the nonmagnetic sites and reduce the magnetic 
moments at magnetic sites, which is an opposite effect to 
a suppression of PD by inducing magnetic moments at 
the nonmagnetic sites. Such study beyond the mean-field 
approximation remains an interesting problem for future 
study. 

To summarize, we have studied the ground state of the 
periodic Anderson model on a triangular lattice by the 
mean-field approximation up to the twelve-site unit cell. 
We have discovered partially disordered states at partic- 
ular commensurate fillings, which are different from the 
previously- found half-filling one. In the new partially- 
disordered states, the nonmagnetic states have a spatially 
extended nature and form an extensive network with 
magnetic sites. This is distinct from the situation at half 
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With further increasing U or decreasing V, the system ex- 
hibits a phase transition with a magnetic moment induced 
at the nonmagnetic site, and finally, changes into the 120° 
antiferromagnetic metallic state continuously. 



